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The  implementation  of  the  BKW  computer  code  on  the  CSIRO  Cyber  76 
computer  network  is  described.  A  user  interface  is  developed  so  that 
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and  so  that  results  can  be  obtained  reliably  and  quickly  from  data  already 
stored  in  the  network. 
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DETONATION  PARAMETERS  FOR  AUSTRALIAN  H I GH - EXPLOS I VES  COMPOSITION 


1 .  INTRODUCTION 


Hiqh  explosives  are  a  fundamental  ingredient  in  most  weapons  and 
devices  used  by  the  Australian  Defence  Force.  Some  explosives  are 
manufactured  in  Australia  while  others  are  imported.  Hiqh  explosives  may  he  a 
sinqle  chemical  compound  (such  as  TNT)  or  a  mixture  (such  as  Composition  Ft). 
Even  for  Composition  B,  the  explosive  components  are  mixed  in  different 
proportions  in  Australia  from  those  used  overseas. 


The  choice  of  explosive  depends  on  its  sensitivity  durinq 
manufacture,  storaqe  and  transport,  its  performance  and  of  course  its  cost. 
Different  environments  (such  as  underwater  in  the  warhead  of  a  torpedo  or  in 
an  air-burst  artillery  shell)  dictate  different  explosive  compositions. 

Systems  analysts  require  data  on  explosive  parameters  so  that  wpsnon 
performance  can  be  predicted.  If  necessary,  the  weapon  design  can  be  modified 
to  optimise  the  terminal  effectiveness. 


The  explosive  parameters  of  interest  in  such  an  analysis  are 
detonation  velocity,  pressure,  temperature  and  enerqy  output.  Their 
experimental  measurement  is  usually  slow  and  potentially  hazardous.  Moreover, 
the  accuracy  to  which  these  measurements  are  made  could  result  in  erroneous 
conclusions  particularly  when  similar  compositions  are  being  compared. 

Theoretical  estimates  of  detonation  parameters  are  based  on  suitable 
high-pressure  equations -of -state  for  the  qaseous  detonation  products.  The 
three  most  commonly  used  are  BICW  (1),  JWL  [21  and  JCZ  [31. 


Of  these,  BKW  has  been  employed  extensively  both  within  Australia 
and  overseas  for  the  last  fifteen  years.  Not  only  has  it  generated  detonation 
parameters  so  that  explosive  performance  can  he  assessed,  but  has  provided 
input  data  for  large  hydrodynamic  codes  used  in  weapon  studies.  Two  codes  in 
operation  in  Australia  are  the  one-  and  two-dimensional  Laaranoian  programs, 
SIN,  [5,61  and  2DL  [71. 


The  first  BKW  proqram  [4]  available  in  Australia  was  written  for  a 
batch -oriented  process  on  a  CDC  3600  computer  althouah  the  earlier  versions 
were  compiled  in  a  machine  languaqe  for  an  IBM  7030  machine.  Followinq  its 
modification  to  permit  runninq  on  the  Cyber  76  computer  in  the  CSIRO  computinq 
network  (CSIRONET),  an  interactive  proqram  USERBKW  was  developed.  By  this 
approach,  the  input  data  deck  is  compiled  directly  from  a  series  of  questions 
and  from  data  already  available.  Results  can  then  be  obtained  reliably  and  at 
minimum  cost. 


This  Report  describes  USERBKW  and  presents  a  table  of  detonation 
parameters  of  high  explosive  compositions  of  current  interest  in  Australia. 


2.  THE  BKW  PROGRAM 


The  BKW  Code  (41  employs  the  Becker-Kistiakowsy-Wi Ison  (BKW) 
equation  of  state  [8]  to  calculate  the  detonation  properties  of  explosives. 

It  is  also  possible  from  these  calculations  to  obtain  the  Huqoniot,  the 
isentrope  throuqh  the  C-J  point  and  the  coefficients  to  fits  of  pressure, 
energy  and  temperature  alonq  the  isentrope.  The  fits  of  these  results  are 
used  in  reactive  hydrodynamic  codes  to  calculate  the  appropriate  equations  of 
state. 


The  BKW  equation  of  state  is  normally  used  in  the  empirical  form: 


PV 

RT 


+  Y 


0Y 

e 


where  y  =  k[v(T  +  0)  ] 
K  =  <  I .  x.  k . 


k.  =  10.46  v. 

i  l 


P  being  the  pressure 

V  the  molar  gas  volume 

R  the  gas  constant 

T  the  absolute  temperature 

x^  the  mole  fraction  of  component  i 

k.  is  the  covolume 

v.  is  the  effective  volume  (in  A®^)  occupied 

by  a  rotating  molecule. 


The  constants  a,  0,  <  and  8  are  determined  experimentally,  and  are  discussed 
more  fully  elsewhere  (8). 


Steady-state  time-independent  behaviour  is  assumed  in  the  use  of 
this  code.  Detonation  products  are  considered  to  be  in  instantaneous  chemical 
equilibrium.  In  practice,  however,  the  detonation  process  is  a  complex  time- 
dependent  one  and  appropriate  caution  in  the  use  of  the  results  is  naturally 
required. 


2 


The  code  itself  contains  a  numher  of  subroutines  that  operate  under 
the  control  of  the  main  proqram.  This  main  routine  handles  only  the  input  and 
output  of  data  besides  the  control  of  the  code.  All  computations  are 
performed  in  the  subroutines.  Ibe  proqram  requires  the  calculation  to  proceed 
iteratively  by  minimisinq  the  free  enerqv. 


3.  SOURCES  OF  DATA 


For  the  explosive,  the  necessary  input  data  are  elemental 
composition,  heat  of  formation,  density  and  formula  weiqht.  For  an  explosive 
mixture,  the  fraction  by  weiqht,  the  elemental  composition  and  heat  of 
formation  of  each  component  are  required. 


The  heats  of  formation  for  some  substances  are  listed  in  various 
publications  [8,9,10,11,12]. 


For  other  compounds,  reliable  estimates  may  be  obtained  by  assuminq 
that  the  heat  of  formation  is  independent  of  the  path  by  which  the  reaction 
proceeds.  Under  this  approach,  the  elements  in  their  standard  states  are 
assumed  to  pass  throuqh  the  monatomic  qaseous  elements  to  the  product 
compound.  Bond  enerqies  are  reasonably  well-known,  and  BKW  is  not  very 
sensitive  to  the  heat  of  formation  used  to  describe  the  explosive. 


The  parameters  a,  0,  ic  and  0  have  been  carefully  calibrated.  In 
qeneral,  the  set  of  parameters  which  best  fits  RDX  (o  =  0.5,  0  =  0.16, 

<  =  10.91,  0  =  400)  are  used  [8].  If  the  explosive  composition  is  deficient 
in  oxyqen  and  produces  siqnificant  concentrations  of  solid  carbon,  the  TNT 
parameters  are  more  appropriate  (a  =  0. 5,  0  =  0.09585,  k  =  12.685  and 
9  =  400). 


The  covolume  for  each  species  is  qiven  by  10.46  v^  where  v^  is  in 
A0"3.  It  represents  the  spherical  volume  of  the  rotatinq  molecule  where  its 
radius  is  the  maximum  dimension  of  the  molecule  as  measured  from  the  centre  of 
mass.  Covolumes  for  the  usual  product  species  are  listed  [41,  [81  or  may  be 
calculated  from  the  Van  der  Waals  radii,  values  of  which  are  tabulated  for 
some  elements. 


The  BKW  code  requires  the  entropy  of  the  product  species  to  be 
expressed  as  a  polynomial  function  of  temperature.  While  these  data  are 
available  for  most  common  species  [8],  information  for  other  products  is 
available  throuqh  tables  [111,  [121  or  throuqh  the  TDF  proqram  [131. 


As  BKW  uses  basically  iterative  techniques,  some  parameters  have  to 
be  selected  initially  so  that  converqence  durinq  the  proqram  takes  place.  For 
example,  initial  estimates  of  pressure  and  temperature  must  be  of  the  riqht 
order.  BKW  is  written  to  ensure  realistic  initial  values  are  selected  for 
CHNO-type  (RDX,  TNT)  explosives.  If  an  appreciably  different  explosive  is 
used,  new  values  of  these  parameters  must  be  selected.  Typical  values  of 
these  modified  parameters  are  qiven  elsewhere  [14], 


For  example,  if  there  is  an  error  in  the  qas  volume  VR0S(1)  (at 
parameter  2)  should  be  chanqed.  If  the  Huqoniot  temperature  is  incorrect, 
HUGBOS  (at  parameter  5)  should  be  amended.  APGCJ  (parameters  11)  or  BPGCJ 
(parameter  12)  should  be  altered  if  there  are  too  many  iterations  in 
calculating  the  detonation  pressure. 

It  may  be  noted  that  aluminium  containing  explosives  tend  to  qive  a 
much  higher  temperature  than  non-aluminised  explosives 


4.  THE  USER  INTERFACE 


The  oriqinal  BKW  proqram  (41  reauired  as  input  a  deck  of  punched 
cards  in  a  complicated  format  specifying  the  density,  the  components  of  the 
explosive  mixture,  the  detonation  products  and  a  larqe  number  of  physical  and 
numerical  constants  related  to  them.  So  complicated  is  the  input  that  is 
requires  a  knowledge  of  the  BKW  proqram,  the  editor  proqram,  the  FORTRAN 
languaqe  and  the  reference  books  for  obtaining  the  numerical  constants.  Data 
preparation  is  both  tedious  and  error  prone,  since  the  input  is  a  precise 
sequence  of  spaces  and  numbers.  There  is  thus  a  clear  need  for  a  user 
interface  packaqe,  such  as  USERBKW,  to  allow  users  with  a  minimum  of  technical 
knowledqe  to  prepare  data  for  BKW  in  a  flexible,  reliable  and  efficient 
manner. 

USERBKW  is  written  as  a  FORTRAN  proqram  for  the  CSIRONET  Cyber  7600, 
although  it  can  be  readily  modified  for  other  computers.  It  runs  under  the 
CIO  interactive  sub-system,  ie  it  asks  questions  of,  and  reads  answers  from, 
the  user  sitting  at  a  remote  terminal.  USERBKW  also  maintains  three  database 
files,  namely  RKWCMP ,  BKWfiAS  and  RKWSOI,.  RtWCMP  contains  information  on 
standard  components  of  explosives,  while  RKWGAS  and  RKWSOL  contain  information 
on  gaseous  and  solid  detonation  products.  These  files  can  be  updated  from 
USERBKW,  and  their  internal  formats  are  given  in  Appendix  T.  Appendix  II 
provides  the  listing. 

The  pies'  ions  ask  e  j  by  " -FRBKW  are  •;«  ra  i  ght  forward  to  answer, 
rxiestions  requiring  a  "yes  .  r  n  "  *  vne  answer  have  the  characters  (Y/N)  at  the 
end  of  the  question.  "Pie  1<"  ii1'  in-wr  '*’iat  [S,  the  effective  answer  if 
nothing  but  a  line  *  ,  "r.  When  asked  to  specify  a  component 

or  a  product ,  tn>.  answer  re  f..  i  ,  i  name  .  t  up  to  six  characters.  For 

example,  TNT  is  a  '•.■.mpon.a  •  ,  u .  i  m>-,,us  product  and  SOliC  is  a  solid 
product.  Typing  a  igest  ion  mar'  ,  r.  u,  we  r  •  • ,  such  a  question  evokes  a  list  of 
standarl  components  ,r  ia  •  di  -t  -  ■,  •  ,  ,  e.  The  urogram  will  automatical  lv  look 
up  the  appropriite  t  -a  mu  la  in  1  ia.v- .  :  ••  ais*ants  if  on  file;  otherwise  the 

user  is  require)  t.  ,  supply  all  suet  !af  i. 

In  answer  to  t  r  i .  mention  "what  n  His  formula  of  ....",  a  chemical 
formula  must  lie  given.  The  formula  must  tie  a  sequence  of  element  svmhols, 
each  of  which  is  followed  by  a  number  or  spaces  or  both.  Brackets  may  be 
used.  For  example,  the  formula  for  hv  irons  emptier  sulphate  is  CU  S  04 
(H20)5.  Flement  symbols  must  not  be  run  together  as  in  CO  which  is  cobalt  not 
carbon  dioxide.  Formulae  should  not  be  confused  with  names.  For  example,  the 


gaseous  product  carbon  dioxide  has  formula  C  02  but  has  name  C02.  If  common 
element  symbols  are  used  in  the  formula,  the  proqram  will  automatically  work 
out  the  formula  weight. 

Some  questions,  (such  as  askinq  for  the  percentage  composition) 
require  a  numerical  answer.  Others  may  represent  a  number  of  alternatives  and 
require  the  user  to  type  a  sinqle  letter  (and  line  feed)  to  indicate  his 
choice.  Examples  of  the  latter  are  the  menus  which  allow  the  user  to  inspect 
and  edit  data  at  will.  There  are  three  such  menus,  the  main  menu,  the 
component  menu  and  the  product  menu. 

The  first  question  that  USERBKW  asks  is  whether  the  user  wants  to 
read  back  in  a  previously  generated  deck.  The  answer  to  this,  under  normal 
circumstances,  is  "no".  The  next  question  asks  the  user  to  type  a  title  of  up 
to  62  characters.  This  title,  toqether  with  the  date,  will  appear  on  each 
page  of  the  BKW  output.  The  next  question  will  ask  the  user  for  the  loadinq 
density  of  the  explosive  (e.o.  2  g/cm  ). 

After  this  the  user  is  requested  to  specify  a  component  of  the 
explosive.  If  the  specified  component  is  not  on  file,  the  user  will  be 
required  to  type  in  its  formula  and  heat  of  formation.  In  any  case,  the  user 
will  have  to  qive  the  percentage  of  explosive  by  mass.  If,  when  all 
components  have  been  specified,  the  sum  of  the  composition  does  not  equal 
100%,  the  proqram  will  qo  to  the  component  menu  for  the  user  to  edit  the 
components  until  it  does. 

After  the  components  have  been  considered,  the  detonation  produ<-  s 
must  be  specified  and  then  their  amounts  estimated.  Usually,  if  the  exp)  ive 
is  a  C-H-N-0  explosive,  this  is  done  automatically;  otherwise  the  det  on 

products  will  be  explicitly  requested,  and  then  estimates  will  be  aske<. 

There  are  three  methods  of  makinq  estimates,  explicitly,  by  which  the  user 
qives  his  estimates,  hierarchically,  by  which  the  user  qives  priorities  with 
first  priority  beinq  qiven  to  the  product  most  likely  to  be  formed  (e.q. 
priority  1  =  H20,  priority  2  =  C02);  and  by  default  estimates.  If  an  attempt 
at  estimation  is  indequate,  the  user  is  told.  He  can  either  try  aqain,  or 
temporarily  iqnore  the  error  until  more  products  can  be  added  later  throuoh 
the  product  menu. 

The  next  item  to  be  considered  is  the  set  of  parameters  for  the 
Becker-Kistiakowsky-Wilson  equation  of  state.  After  this,  if  all  data  are 
present  and  correct,  the  user  may  have  the  proqram  generate  an  "obiect  deck", 
which  is  the  main  output  of  the  package.  The  obiect  deck  contains  a  complete 
job  for  runninq  under  the  ED  subsystem.  This  iob  "chains"  to  the  BKW  proqram 
and  submits  all  appropriate  data.  If  the  user  answers  "yes",  to  the  question 
"do  you  want  to  generate  a  deck  at  this  staqe?"  the  object  deck  is  qenerated 
as  a  file  after  the  user  qives  a  file  name  and  a  user  identification  for  it. 
Otherwise  the  user  may  inspect  and  edit  the  data  by  menu  operations,  and 
generate  an  object  deck  later. 


User  at 
terminal 


Object 

deck 


BKW 

print-out 


Figure  1.  USERBKW  is  the  "frontend"  of  a  program  "chain''. 


A  maximum  of  12  components  mav  comprise  an  explosive  mixture,  which 
cannot  have  more  than  10  elements.  There  can  he  up  to  20  gaseous  and  5  solid 
detonation  products,  none  of  which  may  contain  an  element  not  present  in  the 
mixture.  If  editing  data  in  the  component  or  product  menus,  changes  such  as 
changing  the  composition,  or  adding  another  detonation  product,  will 
necessitate  re-estimation  of  the  product  quantities.  Similarly,  removal  of  a 
component  will  reguire  reconsideration  of  the  detonation  Droducts. 

The  USERBKW  user  interface,  while  working  closely  with  the  bkw 
program,  was  deliberately  written  as  a  separate  program.  This  is  desirable. 
The  software  is  easier  to  manage  for  the  logically  distinct  tasks  of  data 
preparation  and  numerical  processing,  and  it  is  more  convenient  to  run  as 
different  programs.  USERBKW  takes  fractions  of  a  central  processor  second  to 
execute,  and  uses  a  high  priority,  whereas  BKW,  like  most  "number  crunching" 
programs,  consumes  relatively  large  amounts  of  CPU  time,  and  should  therefore 
be  run  at  the  cheaper  non-prime  priority. 

The  user  will  find  that  the  package  is  simple  to  grasp,  fast  to 
operate,  and  that  it  achieves  the  goal  of  facilitating  reliable  calculations 
through  the  BKW  code. 


5.  RESULTS 


The  input  data  for  a  BKW  run  on  Australian  Composition  B  (RDX  55%, 
TNT  45%)  at  a  density  of  1.65  is  shown  in  Appendix  III.  The  corresoondi no 
output  showing  detonation  parameters,  isentrope  and  Hugoniot  and  nirve-f i ttina 
coefficients  are  given  in  Appendix  IV. 

Detonation  parameters  for  other  explosive  compositions  of  interest 
in  Australia  are  listed  in  Table  I. 
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TABLE  1 


DETONATION  PARAMETERS 


Composition 

Density 

q/cm^ 

Gamma 

Det.  Temp. 

K 

Det.  Vel. 

m/s 

Det.  Pressure 

GPa 

TNT 

1.56 

2.81 

2980 

6680 

1  8.  3 

TNT 

1.64 

2.84 

2920 

6920 

20.4 

RDX/wax  92/8 

1.43 

2.84 

2970 

7320 

20.  2 

RDX/wax  92/8 

1.57 

2.90 

2790 

7810 

24.5 

RDX/wax  92/8 

1.67 

2.94 

2630 

8190 

28.4 

Composition  B 

1.65 

2.  92 

2860 

7760 

25.3 

H6  (Note  1  ) 

1 .80 

3.16 

5300 

7370 

23.5 

Ammonium  Picrate 

1.42 

2.92 

2590 

6780 

16.6 

Ammonium  Picrate 

1.50 

2.96 

2500 

7070 

18.9 

Ammonium  Nitrate 

0.67 

2.01 

5810 

5880 

7.7 

Ammonium  Nitrate 

0.80 

2. 1 1 

5750 

6340 

10.4 

Ammonium  Nitrate 

1.00 

2.26 

5600 

7000 

15.2 

Ammonium  Nitrate 

1.70 

2.76 

4360 

9660 

42.  1 

RDX/wax  91/9 

1.60 

2.92 

2680 

7920 

25.6 

RDX/wax  91/9 

1.55 

2.90 

2760 

7730 

23.7 

PBX  W  106 

1.65 

2.87 

3030 

8200 

28.6 

PBX  W  106 

1.68 

2.89 

2980 

8310 

29.8 

PBX  W  106 

1.715 

2.90 

2920 

8440 

31 . 3 

RDX/TNT  60/40 

1.715 

2.94 

2760 

8050 

28.2 

Amatol  50/50  (Note  2) 

1.00 

2.43 

3950 

6430 

1  2.1 

Amatol  50/50 

1.20 

2.74 

2530 

6230 

12.5 

NQ/TNT  50/50  (TNT) 

1.594 

2.68 

3830 

7860 

26.8 

NQ/TNT  60/40  (TNT) 

1.60 

2.67 

4000 

8080 

28.  5 

NQ/TNT  60/40  ( RDX ) 

1.60 

2.  77 

3860 

8280 

29. 1 

NQ/TNT  50/50  (RDX) 

1.594 

2.80 

3700 

8070 

27.  3 

NQ/TNT  40/60  (RDX) 

1.586 

2.82 

3540 

7840 

25.6 

NQ/TNT  40/60  (TNT) 

1 . 586 

2.72 

3670 

7640 

24.9 

Nitroquanidine 

1.65 

2.88 

2890 

8690 

32.  1 

Nitroquanidine 

1.715 

2.92 

2760 

8950 

35.0 

Amatol  50/50  (TNT) 

1.59 

2.  84 

2250 

7280 

22.0 

Amatol  50/50  (RDX) 

1 . 59 

3.02 

2060 

7700 

23.5 

Amatol  80/20 

1.46 

3.02 

1  740 

7490 

20.4 

Torpex 

1  .72 

2.90 

5271 

7150 

22.5 

Hexal  17 

1 . 74 

3.01 

4950 

7710 

25.8 

Nitroquanidine 

1.62 

2.87 

2960 

8570 

30.7 

Astrolite 

1  .  30 

2.99 

1660 

7780 

1  9.  7 

Astrolite 

1.41 

3.04 

1  500 

8300 

24.0 

RDX/TNT/AN  62/10/28 

1.625 

2.96 

2390 

8100 

27.0 

RDX/AN  48.1/51.9 

1.68 

3.10 

1878 

8550 

30.0 

PETN/AN  66.5/33.5 

1.60 

2.  94 

2440 

8000 

25.9 

•  « 


•  « 


•  I 


•  < 


•  i 
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TABLE  1 


(Continued ) 


RDX/TNT/AN 

48.5/10/40.5 

1.670 

RDX/TNT/AN 

25/1 0/65 

1.678 

Tetryl 

1.50 

Metabel  (Note  3) 

1 .608 

Note  1  H6 :  RDX  43%,  TNT 
Note  2  Amatol:  ammonium 
Note  3  Metabel:  PETN  70% 


3.03 

21  20 

8360 

29.0 

3.09 

1730 

8530 

29.8 

2.83 

3170 

6970 

19.0 

2.75 

3010 

7330 

23.0 

30%,  A1 

22%,  NC  0.6%, 

Wax  4.4%. 

nitrate, 

TNT. 

,  TNT  4% 

,  DNT  20%,  NC 

3%  Wax  3%. 

General  Note: 


(RDX)  or  (TNT)  denotes  that  the  calculations  were  carried  out 
with  RDX  or  TNT  parameters. 


6.  CONCLUSION 


BKW  has  proved  a  most  successful  tool  in  predictinq  detonation 
parameters.  It  is  fast,  efficient,  reliable  and  cheap  to  use.  USERBKW  has 
simplified  the  operation  of  the  BKW  code. 
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APPENDIX  I 


MAINTENANCE  OF  USERBKW 


The  FORTRAN  EXTENDED  source  listinq  for  USERBKW,  as  qiv- 
Appendix  II,  is  well  documented  with  comments,  and  the  maintenance 
should  be  familiar  with  these  before  makinq  any  changes.  The  pro- 
under  the  CSIRONET/SCOPE  2.1  operatinq  system.  As  with  most  1-0 
proqrams,  it  relies  on  its  environment.  Subroutines  PAGE  and  PSS' 
instance,  depend  on  the  terminals  in  use,  while  other  parts  of  th- 
make  calls  to  the  operatinq  system  for  file  cataloquinq  operation 

USERBKW  does  a  lot  of  error  checkinq  for  the  benefit  of 
and  may  often  issue  messaqes  in  the  form  of  "**ERROR** . . . " .  Such 
for  the  information  of  the  user.  Very  rarely,  thouqh,  a  diaqnost 
the  form  "=DIAG=..."  may  be  issued.  Such  an  occurrence  may  indie 
presence  of  a  buq,  and  should  he  drawn  to  the  attention  of  the  ma 
programmer. 


The  three  database  files  BKWCMP  (for  components),  RKWGA 
products)  and  BKWSOL  (solid  products)  should  be  maintained  by  the 
within  the  proqram  (e.q.  by  using  the  file  option  in  the  componen 
menus) .  Entries  in  these  files  may  be  created,  deleted  or  update 
entries  may  be  protected  against  accidental  deletion  by  the  use  o 
from  within  USERBKW.  All  users  are  encouraged  to  create  entries  > 
temporarily,  where  this  will  circumvent  tedious  retypinq  of  the  m- 

BKWCMP  is  loqical  file  TAP El ,  and  contains  the  followin' 
about  components:  name,  comment,  formula,  protection,  heat  of  fo 
formula  weight.  Each  entry  is  a  single  line  in  the  format  (1H0, 
2G1 2.6) . 


BKWGAS  is  logical  file  TAPE2  and  contains  the  followina 
on  gaseous  products:  name,  comment,  formula,  protection,  heat  of 
entropy  constants  A,  B,  C,  D,  E  and  IC,  and  the  covolume.  Each  e 
the  format  (1H0,  A6 ,  4A10,  6A10,  A1 ,  G1 2.6/1 X,  7G14.8),  which  tak 
lines. 


BKWSOL  is  TAPE3  with  information  on  solid  products  as  f 
name,  comment,  formula,  protection,  heat  of  formation,  formula  we 
constants  A,  B,  C,  D,  E  and  IC,  covolume,  specific  volume  and  the 
Fickett  solid  equation  of  state  parameters  Ag,  Bg,  C^,  ng,  Eg,  A1 
Cj.  Each  entry  is  in  the  format  1H0,  A6 ,  4A10,  6A10,  A1 ,  2G12.6/ 
9G14.8/1X,  9G14.8)  and  uses  three  consecutive  lines. 

For  each  of  these  three  files  a  file  entry  is  deemed  tc 
there  is  the  character  zero  (0)  in  the  first  column  of  the  line, 
function  SEEK).  A  file  usually  ends  with  an  end-of-partition  mar 
by  a  comment  on  the  last  chanqe  made  to  it.  (See  subroutine  SFAI 
information  such  as  the  internal  formats  of  the  files  is  irreleva 
user,  and  is  only  provided  here  for  the  benefit  of  the  maintenanc 
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